Quark Imaging in the Proton 
Via Quantum Phase-Space Distributions 

A.V. Belitsky0 Xiangdong Ji|l| and Feng Yuar0 

Department of Physics, 
University of Maryland 
College Park, Maryland 20742 

Institute for Nuclear Theory, 
University of Washington 
Box 351550, Seattle, WA, 98195 
(Dated: February 1, 2008) 

Abstract 

We develop the concept of quantum phase-space (Wigner) distributions for quarks and gluons in 
the proton. To appreciate their physical content, we analyze the contraints from special relativity 
on the interpretation of elastic form factors, and examine the physics of the Feynman parton 
distributions in the proton's rest frame. We relate the quark Wigner functions to the transverse- 
momentum dependent parton distributions and generalized parton distributions, emphasizing the 
physical role of the skewness parameter. We show that the Wigner functions allow to visualize 
quantum quarks and gluons using the language of the classical phase space. We present two 
examples of the quark Wigner distributions and point out some model-independent features. 
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I. INTRODUCTION 



In exploring the microscopic structure of matter, there are two frequently-used ap- 
proaches. First, the spatial distribution of matter (or charge) in a system can be probed 
through elastic scattering of electrons, or photons, or neutrons, etc. The physical quantity 
that one measures is the elastic form (structure) factors which depend on three-momentum 
transfer to the system. The Fourier transformation of the form factors provides direct 
information on the spatial distributions. The well-known examples include the study of 
charge distribution in an atom and the atomic structure of a crystal. The second approach 
is designed to measure the population of the constituents as a function of momentum, or 
the momentum distribution, through knock-out scattering. Here the well-known examples 
include the nucleon distributions in nuclei measured through quasi-elastic electron scatter- 
ing, and the distribution of atoms in a quantum liquid probed through neutron scattering. 
The scattering cross section sometimes depends on the reaction dynamics which must be 
understood before the momentum distribution can be extracted. 

Both approaches are complementary, but bear similar drawbacks. The form factor mea- 
surements do not yield any information about the underlying dynamics of the system such 
as the speed of the constituents, whereas the momentum distribution does not give any in- 
formation on the spatial location of the constituents. More complete information about the 
microscopic structure lies in the correlation between the momentum and coordinate spaces, 
i.e., to know where a particle is located and, at the same time, with what velocity it travels. 
This information is certainly attainable for a classical system for which one can define and 
study the phase-space distribution of the constituents. For a quantum mechanical particle, 
however, the notion of a phase-space distribution seems less useful because of the uncertainty 
principle. Nonetheless, the first phase-space distribution in quantum mechanics was intro- 
duced by Wigner in 1932 and many similar distributions have been studied thereafter. 
These distributions have been used for various purposes in very diverse areas such heavy- 
ion collisions, quantum molecular dynamics, signal analysis, quantum information, optics, 
image processing, non-linear dynamics, etc. 2| In certain cases, the Wigner distributions can 
even be measured directly in experiments J 0, 1^ , providing much information about the 
dynamics of a system. 

The main interest of this paper is about the internal structure of the proton (or neutron), 
for which the underlying fundamental theory is quantum chromodynamics (QCD). With 
some changes to accommodate the relativistic nature of the problem, both experimental 
approaches alluded to above have been successfully used to unravel its quark and gluon 
structure: The elastic form factors of the proton have been measured since the 1950s and, 
at low-momentum transfer (< nucleon mass Mat) where the nucleon recoil effects are small, 
the three-dimensional (3D) Fourier transformation of these form factors can be interpreted 
as spatial charge and current distributions of quarks 0. Feynman parton distributions, 
measurable in high-energy inelastic scattering such as deep-inelastic scattering (DIS) and 
Drell-Yan process, have a simple interpretation as the momentum distributions of the quarks 
and gluons in the infinite momentum frame (IMF) 0]. However, the notion of correlated 
position and momentum distributions of quarks and gluons has not been systematically 
investigated in the field, although it is clear that the physics of a phase-space distribution 
must be very rich. 

In this paper, we explore to what extent one can construct physically-interesting and 
experimentally-measurable phase-space distributions in QCD, and what information it con- 
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tains about the QCDparton dynamics. [A brief account of some of the results can be found 
in Ref. j^, see also p].] To facilitate the construction, we examine the uncertainty in the 
traditional interpretation of electromagnetic form factors due to relativity, and analyze the 
physical content of the Feynman parton distributions in the rest frame of the proton. We 
then introduce the phase-space Wigner distributions for the quarks and gluons in the pro- 
ton, which contain most general one-body information of partons, corresponding to the full 
one-body density matrix in technical terms. After integrating over the spatial coordinates, 
one recovers the familiar transverse- momentum dependent parton distributions On the 



other hand, some reduced version of the distributions is related, through a specific Fourier 
transformation, to the generalized parton distributions (GPDs) which have been studied 
extensively in the literature in recent years [l^, EE Ej; EE EE E3 ■ Roughly speaking, 
a GPD is a one-body matrix element which combines the kinematics of both elastic form 
factors and Feynman parton distributions, and is measurable in hard exclusive processes. 
Therefore, the notion of phase-space distribution provides a new 3D interpretation of the 
GPDs in the rest frame of the proton. There are other interpretations of the GPD in the 
literature which are made in IMF and impact parameter space [iEEi,l2El2il23. 



The presentation of the paper is as follows. In Section II, we examine the constraints 
on the physical interpretation of the form factors and parton distributions from relativistic 
effects, anticipating their extension to a full phase-space distribution. In Section III, we first 
briefly summarize the main features of a quantum mechanical Wigner distribution, then 
introduce the quantum phase-space distributions for the quarks and gluons in a rest-frame 
proton. In Section IV, we exhibit the spatial 3D images of quarks generated from slicing the 
quantum phase-space (Wigner) distributions at different Feynman momentum and comment 
on their general features. Section V contains the summary and conclusion. 



II. RELATIVITY CONSTRAINT ON INTERPRETATION OF FORM FACTORS 
AND PARTON DISTRIBUTIONS 

In the literature, the quantum phase-space distributions have been mostly applied to non- 
relativistic systems. For the proton, however, relativity is essential. In measuring the elastic 
form factors of the proton, the momentum transfer to the system can easily exceed the rest 
mass, resulting a large recoil and Lorentz contraction. The quarks and gluons inside the 
proton follow relativistic dynamics. Moreover, when a quark is struck in a DIS experiment, 
it travels along the light-cone: the trajectory of an extreme-relativistic particle. Therefore, 
to develop a phase-space distribution of the proton, we must examine to what extent the 
notion actually makes sense for relativistic systems. 

In the first subsection, we examine the textbook interpretation of the electromagnetic 
form factors of the proton, reminding the reader that there are intrinsic ambiguities in the 
interpretation. We emphasize, however, that different ways of the interpreting the form 
factors can be regarded as different choices of schemes. When used consistently, one scheme 
is in principle as good as any other. The degree of scheme dependence depends on l/(Mi?), 
where M is the mass and R is some kind of radius, which is 1/4 for the proton. 

In the second subsection, we consider the Feynman parton distributions, most-commonly 
interpreted as the momentum densities in IMF. Since the notion of a phase-space distribu- 
tion is meant for a proton in its rest frame, and since the distribution should be reduced to 
the Feynman distribution after integrating out the spatial coordinates, we are compelled to 
examine the physics of latter in the static system of coordinates. In particular, we need to 
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understand the meaning of Feynman momentum x in that context. This can be achieved by 
introducing the so-called spectral function — the correlated momentum and energy distribu- 
tion of the constituents — familiar in non-relativistic many-body physics. In the process, we 
find that the separation between particle and antiparticle, familiar in the IMF, disappears: 
One can only keep track of the creation and annihilation of fermion quantum numbers, such 
as the electric charge. 



A. The Proton Form Factors and Scheme-Dependent Charge Distributions 

The electromagnetic form factors are among the first measured and mostly studied ob- 
servables of the nucleon. They are defined as the matrix elements of the electromagnetic 
current between the nucleon states of different four-momenta. Because the nucleon is a spin 
one-half particle, the matrix element defines two form factors, 

(P2|j,(0)bi) = f/fe) [F,{q')l, + f/(pi) , (1) 

where Fi and F2 are the well-known Dirac and Pauli form factors, respectively, depending on 
the momentum transfer q = P2 — Pi, and U{p) is nucleon spinor normalized as U{j>)U{p) = 
2Mn. 

Since the beginning, it has been known that the physical interpretation of the nucleon 
form factors is complicated by relativistic effects . Consider a system of size R and mass 



M. In relativistic quantum theory, the system cannot be localized to a precision better than 
its Compton wavelength 1/M. Any attempt to do this with an external potential will result 
in creation of particle-antiparticle pairs. As a consequence, the static size of the system 
cannot be defined to a precision better than 1/M. If R ^ '^/M, which is the case for all 
non-relativistic systems, the above is not a significant constraint. One can probe the internal 
structure of the system with a wavelength comparable to or even much smaller than 

R, but still large enough compared to 1/M so that the probe does not induce an appreciable 
recoil. A familiar example is the hydrogen atom for which RMh ~ Mh /{meaem) ~ 10^, 
and the form factor can be measured through electron scattering with momentum transfer 

When the probing wavelength is comparable to 1/M, the form factors are no longer 
determined by the internal structure alone. They contain also the dynamical effects of 
Lorentz boosts because the initial and final protons have different momenta. In relativistic 
quantum theory, the boost operators involve nontrivial dynamical effects which result in 
the nucleon wave function being different in different frame (in the usual instant form of 
quantization). Therefore in the region |g| ~ M, the physical interpretation of the form 
factors is complicated because of the entanglement of the internal and the center-of-mass 
motions in relativistic dynamics. In the limit |g| ^ M, the former factors depend almost 
entirely on the physical mechanism producing the overall change of the proton momentum. 
The structural effect involved is a very small part of the nucleon wave function (usually the 
minimal Fock component only). 

For the nucleon, MjsiRj^i ~ 4. Although much less certain than in the case of the hydrogen 
atom, it seems still sensible to have a rest-frame picture in terms of the electromagnetic form 
factors, so long as one keeps in mind that equally justified definitions of the nucleon sizes 
can differ by ~ 1/ Mn{RnMn)- For example, the traditional definition of the proton charge 
radius in terms of the slope of the Sachs form factor ^^(g^) is 0.86 fm. On the other hand, 
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if one uses the slope of the Dirac form factor Fi to define the charge radius, one gets 0.79 
fm, about 10% smaller. The form factors at |g| > Mn ~ 1 GeV cannot be interpreted as 
information about the internal structure alone. 

To further clarify the uncertainty involved in the interpretation of the electromagnetic 
form factors, let us review the textbook explanation offered originally by Sachs 

To establish the notion of a static (charge) distribution, one needs to create a wave-packet 
representing a proton localized at R 

\R) = J ^s^'^'^^'^iP^lP)^ (2) 

where the plane-wave state \p) is normalized in a relativistic-covariant manner {p2\pi) = 
2Ep_^{27iYS^^\pi — P2), and ^{p) is the momentum space profile normalized as 
2 / d^pEpl^lpjl"^ = (27r)'^. The wave-packet is not an eigenstate of the free Hamiltonian. 
Therefore, as time progresses, the wave packet will spread. The characteristic spreading 
time is proportional to {Mp^/fp') which is long for a non-relativistic system. But for a 
relativistic particle, the spread could happen much faster compared to the characteristic 
time-scale of a weakly-interacting probe. 

Having localized the wave-packet at i? = 0, we can calculate, for example, the charge 
distribution in the wave-packet 

p(f) = (i? = 0|jo(r)|i? = 0), (3) 

where r measures the relative distance to the center i? = 0. Taking its Fourier transform, 
one gets 



F{q) = J d^fe''^-''p{r 
d^p 



^* {p + q/2) ^{p- q/2) {p + q/2\ jo(0) \p- q/2) , (4) 



where we have changed the momentum integration variables, with p representing the average 
momentum of the initial and final protons. It is important to point out that the resolution 
momentum q is now linked to the difference in the initial and final state momenta. In non- 
relativistic quantum systems, because of the large masses, the momentum transfer causes 
little change in velocity, and hence the initial and final states have practically the same 
intrinsic wave functions. In relativistic systems, this is the origin of the difficulty in inter- 
preting the form-factor: we do not have a matrix element involving the same intrinsic proton 
state. 

To remove the effects of the wave packet, the necessary condition on is that the 
coordinate-space size of the wave-packet must be much smaller than the system size \5r\ -C 
Rn- Furthermore, the probing wave length, or resolution scale, must also be large compared 
with the size of the wave-packet 5r~l/p^l/g. Then one can ignore g'-dependence in 
so that ^ {p± \q^ ^ ^{p) 

m = J -0^,\m\'{p+m\Mo)\p-q/2) ■ (5) 

On the other hand, to be insensitive to the anti-particle degrees of freedom, the size of the 
wave packet must be larger than the proton Compton wave length \6f\ ^ l/Mp^. In the 
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momentum space this corresponds to restrictions on momenta allowed in the wave packet 
|p| -C Mat. Therefore the combined constraint on the wave packet profile is 1/Rn -C 
|g| -C IpI -C Mat. The extreme limit of the last inequality yields a wave packet with a 
zero-momentum nucleon 

l'^(p)P = |^5(^nPl, (6) 

which gives 

2M;vF(g-) = (g72|jo(0)|-g72) • (7) 



This is the matrix element of the charge density in the Breit frame, and is 2MnG E{t)'W2'Wi 
where 

GE{t)=F^{t) + ^F,{t) (8) 

is the Sachs electric form factor {t = —(f) and the Weyl spinors involved are normalized 
conventionally by w*w = 1. Hence, we arrive at the textbook interpretation of G^; as a 
Fourier transformation of the proton charge distribution. 

Likewise, the Sachs magnetic form factor Guit) = -^i(^) + ^2(t) is obtained from the 
Breit frame matrix element of the electric current 

m\ J(0) hq/2) = X g1GM(t) , (9) 

where the three- vector of spin is s*= W2^(twi. 

It must be pointed out that the charge and magnetization distributions thus defined 
contains the Lorentz contraction effects along the photon direction q when (f ^ 4M^, 
which make the proton look like a pancake. Various prescriptions exist in the litera- 
ture which have been proposed to remove the relativity effects and extract the "intrinsic" 
charge/magnetization distributions from the experimental data 0, HME^- However, it is 
difficult to accomplish it in a model-independent way. 

Since relativity makes the interpretation of the electromagnetic form factors non-unique, 
the best one can do is to choose one particular interpretation and work consistently. For 
example, when extracting the proton charge radius from the Lamb shift measurements, one 
shall use the same definition as from the electric form factor. The most frequently-used 
definition is that of Sachs, but other schemes are equally good and the scheme dependence 
disappears in the limit MR oo. This is very much like the renormalization scheme 
dependence of parton densities due to radiative corrections at finite strong coupling constant 
as'. Although the MS scheme is the most popular in the literature, one can use the parton 
densities in any other scheme to correlate physical observables. In this paper, we adopt the 
Sachs interpretation of the form factors. 

Relativistic corrections and Lorentz contraction effects in the transverse dimensions are 
found to disappear in an IMF There the proton has an infinitely large effective mass, 
and hence for physics in the transverse dimensions, we are back to the non-relativistic case. 
In particular, one can localize the proton in the transverse coordinate space with no recoil 
corrections. The Dirac form factor Fi is found to be related to the charge distribution 
in transverse plane, with information along the z-axis integrated. The price one pays for 
eliminating the relativistic effects is the loss of a spatial dimension. 
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B. Parton Distributions As Seen in the Rest Frame of the Proton 



Parton distributions were introduced by Feynman to describe deep-inelastic scattering 
0]. They have the simplest interpretation in the IMF as the densities of partons in the 
longitudinal momentum x. In QCD, the quark distribution is defined through the following 
matrix element, 

g(^) = ^e^'^{pmOh^^i\n)\p) , (10) 

where we have used the standard light-cone notation p"^ = {p'^ ±p^)/\/2, and is a vector 
along the direction of (1, 0, 0, —1) and n ■ p = 1. ip' is a quark field with an associated gauge 
link extending from the position of the quark to infinity along the light cone, and hence is 
gauge-invariant in non-sin gula r gauges. The renormalization scale dependence is implicit. 
In light-cone quantization (271], it is easy to get 

1 f d'k^{p\b{{k+,k^)bxik+,k^)\p) 



2x^t^^J (27r)3 {p\p) 
_ -1 f d'^k± {p\d{{k+,k±)dxik+,k±_)\p) 

where 6^ and d'^ are creation operators of a quark and an anti-quark, respectively, with 
longitudinal momentum k~^ = xp'^ and transverse momentum k±. The interpretation as 
parton densities is then obvious. 

To construct the quantum phase-space distributions for the quarks, we need an inter- 
pretation of the Feynman densities in the rest frame. This is because the IMF involves a 
Lorentz boost along the z-direction which destroys the rotational symmetry of the 3D space. 
However, if one works in the rest frame of the proton, the two quark fields in Eq. (|iup are 
not at the same time. If one Fourier-expands one of the fields in terms of quark creation and 
annihilation operators, the other must be determined from Heisenberg equation of motion. 
The result is that the bi-linear quark operator takes a very complicated expression in terms 
of the creation and annihilation operators in the equal-time quantization. 

The physics of the Feynman quark distribution in the rest frame is made more clear 
through the notion of the spectral function 

s{k) = ^J d'^e''<{p\mi^m\p) ■ (12) 

which is the dispersive part of the single-quark Green's function in the proton. The physical 
meaning of S{k) can be seen from its spectral representation, 

S{k) = ^(27r)^5(^)(p-A;-j9„)(p|^fc|n)7+(n|<Z^(0)|p)/2p+ 

n 

~ 5:(27r)V^)(p-A:-pOI(^l'^.+ b)P (13) 

n 

where \Pk is a Fourier transformation of ^{^): It is the probability of annihilating a quark 
(or creating an antiquark) of /o'ur-momentum k (three-momentum k and the off-shell energy 
E = k^) in the nucleon, leading to an "on-shell" state of energy-momentum Pn = p — k. The 
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quark here is off-shell because if pn and p are both "on-shell" , k"^ ^ rn^ in general. [That 
the partons are off-shell are in fact also true in the IMF calculations.] Of course, in QCD 
\n) is not in the Hilbert space, but the spectral function itself is still a meaningful quantity. 

Since the quarks are ultra-relativistic, \Pk contains both quark and antiquark Fock op- 
erators. One cannot in general separate quark and anti-quark contributions, unlike in the 
non-relativistic systems in which only the particle or antiparticle contribute. In fact, if one 
expands the above expression, one finds pair creations and annihilation terms. However, 
this is also true for the charge density discussed in the previous subsection. Therefore we 
can speak of S{k) as a distribution of vector charges and currents, but not a particle den- 
sity. In nuclear physics where the non-relativistic dynamics dominates, the nucleon spectral 
function in the nucleus is positive definite and can be regarded as a particle density. The 
nuclear spectral function is directly measurable through pick-up and knock-out experiments, 
in which E and k are called the missing energy and missing momentum, respectively (see 
for example [28]). 

It is now easy to see that in the rest frame of the proton, the Feynman quark distribution 

is 

^^''^ = ^/ (0^^^^° + " xMr,)S{k) . (14) 

The X variable is simply a special combination of the off-shell energy k^ and momentum k^. 
The parton distribution is the spectral function of quarks projected along a special direction 
in the four-dimensional energy-momentum space. The quarks with different k^ and k^ can 
have the same x, and moreover, the both x > and x < distributions contain contributions 
from quarks and anti-quarks. 

To summarize, in the proton rest frame, the quarks are naturally off-shell, and hence have 
a distribution in the four-dimensional energy-momentum space. The Feynman distribution 
comes from a reduction of the full distribution along a special direction. 



III. QUANTUM PHASE-SPACE (WIGNER) DISTRIBUTIONS 

In classical physics, a state of a particle is specified by its position r and momentum p. 
In a gas of classical particles, the single-particle properties are described by a phase-space 
distribution f{f,p) representing the density of particles at a phase-space point {f,p). Time 
evolution of the distribution is governed by the Boltzmann equation, or Liouville equation 
if the particles are not interacting. 

In quantum mechanics, position and momentum operators do not commute and hence, 
in principle, one cannot talk about a joint momentum and position distribution of particles. 
Indeed the quantum mechanical wave functions depend on either spatial coordinates or 
momentum, but never both. Nonetheless, Wigner introduced the first quantum phase-space 
distribution just a few years after quantum mechanics was formulated |l|. It is not positive 
definite and hence cannot be regarded as a probability distribution. However, it reduces to 
the positive-definite classical phase-space distribution in /i — limit. The sign oscillation 
in the phase-space is necessary to reproduce quantum interference. The Wigner distribution 
contains the complete single-particle information about a quantum system (equivalent to 
the full single-particle density matrix), and can be used to calculate any single-particle 
observable through classical-type phase-space averages. 

In this section, we first remind the reader some basic features of the quantum phase-space 
(Wigner) functions. We then generalize the concept to the relativistic quarks and gluons in 
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the proton. With the preparation in Section II, the construction is straightforward. How- 
ever, the most general phase-space distribution we define is not measurable at present, and 
hence we proceed to make reductions by integrating out some dependent variables. After 
integrating out the spatial coordinates, we recover the transverse-momentum dependent par- 
ton distributions Upon integrating over the parton transverse momentum, we have the 
reduced Wigner distributions depending on 3-space coordinates and Feynman momentum x, 
which are related to the GPDs by a simple Fourier transformation. Therefore, the reduced 
quantum phase-space distributions are physical observables. 



A. General Aspects of Wigner Distributions 

There is a vast literature on the quantum phase-space distributions, and the Wigner 
distributions in particular. In this subsection, we would like to summarize some of the 
salient features. 

Suppose we have a one- dimensional quantum mechanical system with wave function ip{x), 
the Wigner distribution is defined as 

W{x,p) = J dr]e'P''^*{x - r]/2)tlj{x + r]/2) , (15) 

where we have set ^ = 1. When integrating out the coordinate x, one gets the momentum 
density which is positive definite. When integrating out p, the positive-definite 

coordinate space density follows. For arbitrary p and x, the Wigner distribution 

is not positive definite and does not have a probability interpretation. Nonetheless, for 
calculating the physical observables, one can just take averages over the phase-space as if it 
is a classical distribution 

{0{x,p)) = J dxdpW{x,p)0{x,p) (16) 

where the operators are ordered according to the Weyl association rule. For a single-particle 
system, the Wigner distribution contains everything there is in the quantum wave function. 
For a many-body system, the Wigner distribution can be used to calculate the averages of 
all one-body operators. Sign changes in the phase-space are a hint that it carries non-trivial 
quantum phase information. 

In the classical limit, the Wigner distribution is expected to become classical phase-space 
distribution. For systems which are statistical ensembles, the limit h ^ is often well- 
behaved. For example, for an ensemble of harmonic oscillators at finite temperature, the 
Wigner distribution becomes the classical Boltzmann distribution as ^ 0, see, e.g., j29l |. 
The Wigner distribution for the nth excited state of the one-dimensional harmonic oscillator 
of energy En = hu (n + is js^l 

Wnip, X) = i^e-^^/(ML„ (i^) , (17) 

where H stands for the hamiltonian H{p, x) = p"^/ {2m)+muj'^x^ /2 and L„ is the nth Laguerre 
polynomial. In the quasi-classical limit — vanishing Planck constant and large quantum 
numbers — the oscillator Wigner distribution turns into the generalized distribution resided 
on the classical trajectories = fixed , 

lim Wnip, x) ~ 6(h{p, x) - Ejj . (18) 
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Phase-space averaging with this kernel is equivalent to calculating observables using classical 
equations of motion. This can be easily understood from the semi-classical form of the wave 
function 

ij{x) = C(x)e^^(^)/^ . (19) 

Substituting this into Eq. (fT3j) and expanding S to the first order in h, one gets the quasi- 
classical Wigner distribution 

W{p,x) = \C\'5(p-^^ , (20) 

where the argument of the delta-function describes a family of classical paths. 

The quantum-mechanical Wigner distribution is measurable. The actual measurement 
has been performed for a simplest quantum system — the quantum state of a light mode (a 
pulse of laser light of given frequency) — employing ideas of Vogel and Risken It was 
extracted via the method of homodyne tomograhy by measurement of a marginal ob- 
servable and subsequent reconstruction by the inverse Radon transformation. Recently this 
Wigner distribution has been measured directly by means of the photon counting techniques 
based on a Mach-Zender interferometric scheme [5|. 

Other versions of t he p hase-space distributions are possible. They are all members of the 
so-called Cohen class jsjj , with Husimi and Kirkwood distributions [3^ being its well-known 
representatives. The Husimi distribution is a smeared version of the Wigner distribution 
defined by projection of the wave function on the coherent state (Gaussian wave packet) 



H{p,x) = J dp'dx'W{p',x')W,oh{p' 



p, x' — x) 



which is real and positive-definite. On the other hand, the Kirkwood function is complex. 
All these distributions are expected to reduce to the same phase-space distribution in the 
— s> limit. 



B. Quantum Phase-Space Quark Distributions in the Proton 

In this subsection, we generalize the concept of phase-space distributions to relativistic 
quarks and gluons in the proton. In quantum field theory, the single-particle wave function 
must be replaced by quantum fields, and hence it is natural to introduce the Wigner operator^ 

Wrir,k) = J d\e'^'^'^{f-7]/2)r^{r + 7]/2) , (21) 

where r is the quark phase-space position and k the phase-space four-momentum conjugated 
to the spacetime separation rj. r' is a Dirac matrix defining the types of quark densities 
because the quarks are spin-1/2 relativistic particles. Depending on the choice of we can 
have vector, axial vector, or tensor density. 

Since QCD is a gauge theory, the two quark fields at different spacetime points are not 
automatically gauge-invariant. One can define a gauge- invariant quark field by adding a 
gauge link to the spacetime infinity along a constant four- vector n^, 

^{rj) = exp (^-ig J d\n- A{Xn + 77) j ^'(77) , (22) 
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where we assume the non-singular gauges in which the gauge potential vanish at the space- 
time infinity 3^ 3^, 3^ 3^ 37 1. Clearly, the Wigner operator depends on the choice of n^. 
While theoretically any is possible, experimentally is constrained by the probes. 

We have extended the Wigner distribution to including the time variable. Therefore, 
beside the dependence on the 3-momentum, there is also a dependence on the energy. For 
the bound states in a simple system such as those in a simple harmonic oscillator, the energy 
dependence is a ^-function at the binding energies. For many-body systems, however, the 
energy-dependence is more complicated, as it reflects the distribution of the states after one 
particle is removed from the system. 

For non-relativistic systems for which the center-of-mass is well-defined and fixed, one 
can define the phase-space distributions by taking the expectation value of the above Wigner 
operators in the = state. For the proton for which the recoil effect cannot be neglected, 
the rest-frame state cannot be uniquely defined, as discussed in Section II. Here we follow 
Sachs, defining a rest-frame matrix element as that in the Breit frame, averaging over all 
possible 3-momentum transfers. Therefore, we construct the quantum phase-space quark 
distribution in the proton as. 



Wr{r,k) 



2Mn 
1 



2M, 



N 



q/2 Wr{r,k) - q/2 
"^■''/q/2 >Vr(0,A;) - q/2 



(23) 



where the plane-wave states are normalized relativistically. The most general phase-space 
distribution depends on seven independent variables. 

The only way we know how to probe the single-particle distributions is through high- 
energy processes, in which the light-cone energy k- = (fco - k')/V2 is difficuh to measure, 
where the 2;-axis refers to the momentum direction of a probe. Moreover, the leading ob- 
servables in these processes are associated with the "good" components of the quark (gluon) 
fields in the sense of light-cone quantization which can be selected by r' = 7"*", 7'*'75, or 
(J where 7"*" = (7^+7^) / V^- The direction of the gauge link, n^, is then determined by the 



trajectories of high-energy partons traveling along the light-cone (1, 0, 0, —1) [36|,|37[. There- 
fore, from now on, we restrict ourselves to the reduced Wigner distributions by integrating 
out k~ , 

Wr{r, k) = J ^W^r(r, k) , (24) 

with a light-cone gauge link is now implied. Unfortunately, there is no known experiment 
at present capable of measuring this 6-dimensional distribution which may be called the 
master or mother distribution. 

Further phase-space reductions lead to measurable quantities. Integrating out the trans- 
verse momentum of partons, we obtain a 4-dimensional quantum distribution 



fr{f,k' 



d'^k^ 
(2^ 



Wr{r,k) 



drj 



q/2 n-rmrnrm - m) ■ (25) 



2Mn J (27r)3 J 2tx 

The matrix element under the integrals is what defines the GPDs. More precisely, if one 
replaces k'^ by Feynman variable xp^ (p+ = Eq/\f2^ proton energy Eq = JM^ + cp/A ) 



11 



and rj by A/p"*", the reduced Wigner distribution becomes the Fourier transformation of the 
GPD Frlx, ^, t) 

In the present context, the relation between kinematic variables are ^ = q^/{2Eq) and 
t = —q^. Taking T = v^T"*", the corresponding GPD has the expansion 2 



F,^{x,^,t) = y — e^^^(g72|^(-An/2)£x/27+^(An/2)|-g72 



= H{x,^,t)UimV2^^Ui-m + Eix,^,t)U{m'-^Ui-q/2) , (27) 

where C is the shorthand for the light-cone gauge link. 

The phase-space function (r, x) can be used to construct 3D images of the quarks 
for every selected Feynman momentum x in the rest frame of the proton. These images 
provide the pictures of the proton seen through the Feynman momentum (or "color" or x) 
filters. They also may be regarded as the result of a quantum phase-space tomography of 
the proton. We remind the reader again that the Feynman momentum in the rest-frame 
sense is a special combination of the off-shell energy and momentum along z, namely E + k^. 
Integrating over the z coordinate, the GPDs are set to ^ ~ = 0, and the resulting two- 
dimensional density f'y+{r±,x) is just the impact-parameter-space distribution [l^. Further 
integration over f\L recovers the usual Feynman parton distribution. 

The physical content of the above distribution is further revealed by examining its spin 
structure. Working out the matrix element in Eq. (jTTj) . 

(x, e, t) = [H{x, t) - rE{x, ^, t)] + t[s x g^-i- [H{x, ^, t) + E{x, ^, t)] , (28) 



where r = (f/AMj^. The first term is independent of the proton spin, and is considered as 
the phase-space charge density 

p+(r, x) = J ^e-'n^(a;, t) - rE{x, ^, t)] . (29) 

The second term depends on the proton spin and can be regarded as the third component 
of the phase-space vector current 

m^) = I (^^"^'''^[^"'x ^'^mx,^,t)+Eix,^,t)] . (30) 

The i?-term generates a convection current due to the orbital angular momentum of massless 
quarks and vanishes when all quarks are in the s-orbit. The physics in separating into 
p+ and can be seen from the Dirac matrix 7"*" selected by the high-energy probes, which 
is a combination of time and space components. Because the current distribution has no 
spherical symmetry, the quark charge seen in the infinite momentum frame, p+ -|- j^, is 
deformed in the impact parameter space jHS]- This is the kinematic effect of Lorentz boost. 

Integrating the phase-space charge distribution p^{f,x) over x, one recovers the 
spherically-symmetric charge density in space. On the other hand, if integrating over x 
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in j^{f,x), one obtains the electric current density. In the latter case, if the integral is 
weighted with x, one obtains the mechanical momentum density :8]. 

Finally, when integrating over r in the reduced Wigner distributions in Eq. ()24|) . one 
obtains the transverse-momentum dependent parton distributions. 

There is a lot of interesting physics associated with these distributions which has been 
discussed recently in the literature. For instance, in a transversely polarized proton, the 



quark momentum distribution has an azimuthal angular dependence j39|, |40|, |4l|. The so- 
called Siver's function can produce a novel single-spin asymmetry in deep-inelastic scattering. 
We will not pursue this topic here, except emphasizing that they have the same generating 
functions as the GPDs. 



IV. THREE-DIMENSIONAL IMAGES OF THE QUARKS IN THE PROTON 

Once the GPDs are extracted from experimental data or lattice QCD calculations 
liE S, EE 13 ) the phase-space distributions can be obtained by straightforward Fourier 
transformations. Without a first-hand knowledge on the GPDs at present, we may be able 
to learn some general features of the phase-space distributions form GPD models. 

The GPDs have been parametrized directly to satisfy various constraints [3, 12, 



including 1) the first moments reducing to the measured form factors, 2) the forward limit 
reproducing the Feynman parton distributions, 3) the x-moments satisfying the polynomial- 
ity condition jl^], and 4) the positivity conditions ^S]. In the first subsection, we introduce 
a new parametrization without assuming factorized dependence on the t and other variables. 

The GPDs were first calculated in a realistic nucleon model — the MIT bag model |49|. 
They have also been calculated in the chiral-quark soliton mo del HI |5l|. Recently, there are 



calculations in the quark models as well |52l. I53|. In the second subsection, we will consider 
the Wigner distributions in the bag model. 



A. A GPD Parametrization 

A generalized parton distributions depend on three variables, x, ^, and t. The simplest 
way to satisfy the polynomiality condition is to relate it to a double distribution 0, 4^ 
and the D-term |50] 

H{x,^,t)= j'^^Siy\x,OF[y,^,t^+e{^>\x\)D(^^,t^ , (32) 

The "step" -function kernel in Eq. ()32|) has the form 

s{y\x,o = 0{x > 00 {j^ > y > ^) 

+oi-^ > x)o > y > ^) + > \x\)e > y > ^) . 
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The g-flavor double distribution Fq = F^'^^ + F^'^'^, including both valence and sea, can be 
related to the non-forward quark distribution fq{y,t) through a profile function n{y,z,b), 

F^^\y,z,t) = rAy,t)e{y)7r{\y\,z;K,i), (33) 
^(y, z, t) = {fq{y, t)e{y) - fq{-y, m-y)) 7r(||/|, 6,ea) , (34) 

where at t = the function fq{y,t = 0) reduces to the conventional parton distribution 
functions. The profile function with a single parameter b is assumed to be universal for 
valence- and sea-quark species, and reads j42|] 

.(y z- b) - ^(^ + 3/2) [{l-yr-z^t .35. 



To proceed further, we design a non-factorized ansatz |15|, |2l|, |38| for the functions fg{y, t) 
with intertwined t and y dependence. This is opposed to a factorized form of GPDs with 
completely disentangled dependence of the momentum transfer t and scaling variables (x, £). 



The latter is currently accepted in almost all evaluations of physical observables |17l . 14 71 . 154| . 
Due to a limited kinematical coverage in the t-channel momentum transfer t in experiments, 
theoretical estimates confronted to data are currently insensitive to this feature. Our model 
will be based on the GRV leading order quark distributions [H^ with discarded fiavor asym- 
metry of the sea and it reads 

f:^\y, t) = 1.23%-""-"^(^-^)'''* (1 -1.8^ + 9.5y) (1 - y)'-'' , (36) 
rAy, t) = 0.761y-"" (2y-<(i-^)^'^* - y-z^id-.)*) (1 - 1.8 + 9.5y) (1 - yf'' , 

fuiy,t) = Uy,t)=0.7Qy-'^^-<^'-yy''''il-3.Q^ + 7.8y)il-yA. 

These models naturally reduce to the quark form factors with the dipole parametrization of 
proton and neutron Sachs form factors. The valence ci-quark function has a more complicated 
structure since the corresponding form factor F^ has a node at \t\ ~ 4M^/|2k„ + Kp + 1|: 
it is positive below this value and is negative above it. The Regge intercepts and slope 
parameters are taken as 

ay = 0.52, a; = l.lGeV-^ = 1.0 GeV~% (37) 

a, = 0.85, a; = 0.3GeV-^ 

The valence quarks Regge parameters are numerically close to the ones of p-reggeons, while 
the sea quarks being generated by gluon radiation are analogous to the one of the pomeron. 
The form factor asymptotics at large t is governed by the large-?/ behavior of f{y,t). If the 
latter has the form f{y,t) ~ y-°'-°' (^~y)^^ (^l — then the corresponding form factor is 
F(t — * 00) ~ |j(:|-(^+i)/(p+i). The perturbative QCD asymptotics for valence quarks requires 
p = 1. We use however p = 1/2 for them since this value fits better the form factor at small 
and moderate t. For p = 1 one can get a decent behavior at moderate t with a'^ = 1.6 GeV^. 
We use in our estimates 6vai = &sea = 1- The D-term is parametrized as 

D{z,t)=(^l-^^ {l-z^){d,c!/\z) + ---) , (38) 
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x=0.01 x=0.40 x=0.70 




FIG. 1: The u-quark phase-space charge distribution at different values of the Feynman momentum 
for non-factorizable ansatz of generahzed parton distributions (|^-{6|) . The vertical and horizontal axis 
corresponds to z and |rj_|, respectively, measured in femtometers. The [dashed] contours separate 
regions of positive [darker areas] and negative [lighter areas] densities. Below each contour plot we 
presented the shape of three-dimensional isodensity contours [p = const] . 



with the mass scale = 0.6 GeV^ and the parameter (io computed within the xQSM 
isi 51 1 and on the lattice ji^li^ with the results 

rfXQSM ^ _4 1 ^ ^latt = rfn ^ rf^ ^ _0.5 , (39) 

respectively, where N j is the number of active flavors. In the lattice case, the effect of 
disconnected diagrams was not calculated, however they are known to produce a sizable 
negative contribution 4^1. Once the latter are properly taken into account the lattice result 



might approach the model calculation. For our present estimate we chose an intermediate 
value (io = —1.0. 

According to the previous section, the phase-space charge distribution p+(r, x) is just the 
Fourier transformation of the above GPDs, 

pW. x)=J ^e-^^-^iJ«(x, t) , (40) 
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x=-0.05 



x=0.60 





-1 -0.5 D 0.5 1 



.5 0.5 



-1 -0.5 0.5 1 



-1 -0.5 0.5 1 



FIG. 2: The phase-space charge distribution for the u-quark at negative Feynman momentum 
X = —0.05 and x = —0.4 [two left panels] and d-quark for positive x = 0.4 and x = 0.6 [two right 
panels] . 



where ^ = qz/2Eq, Eg = y M"^ + and t = (f . In the following, we consider the result 
of the quark densities from the above parametrization. 

In Fig.^we show the up-quark charge distributions calculated from Hu{x, ^, t) for various 
values of a; = {0.01,0.4,0.7}. While the intensity of the plots indicates the magnitude of 
the positive distribution, the lighter areas below the ground-zero contours indicate negative 
values. The plots show significant change in the distribution on the longitudinal momentum 
fraction x. The image is rotationally symmetric in the f±-pla.ne. At small x, the distribution 
extends far beyond the nominal nucleon size along the z direction. The physical explanation 
for this is that the position space uncertainty of the quarks is large when x is small, and 
therefore the quarks are de-localized along the longitudinal direction. This de-localization 
reflects a very peculiar part of the nucleon wave function and shows long-range correlations 
as verified in high-energy scattering. In a nucleus, the parton distributions at small x are 
strongly modified because of the spatial overlap between the nucleons. On the other hand, at 
larger x, the momentum along z direction is of order nucleon mass, the quarks are localized to 
within l/M^v. The quantum mechanical nature of the distribution becomes distinct because 
there are significant changes in the sign at different spatial regions. 

It is also interesting to explore the distribution at negative x. We show in Fig. |21 [two left 
panels] the Wigner densities for the up-quark in the proton for x: —0.05 and —0.4. These 
plots show significantly- different pattern than those of the positive x. Finally, the two right 
panels we show the density for the down-quark in the proton. The essential features are 
quite similar to those of the up-quark densities. 



B. The MIT Bag Model 

The MIT bag model was invented more than a quarter of a century ago js^. The model 
was motivated by the color confinement property of QCD. Massless quarks are confined to 
a cavity of radius R, and move freely inside. The quark wave function is ultra-relativistic 
and can be solved from the free Dirac equation with spherical boundary conditions. The 
bag model has been used to calculate many static properties of the nucleon and has had 
many notable successes. The model can also be used to describe the excitation spectrum of 
hadrons :56] . The electromagnetic form factors js^ and parton distributions have also been 
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calculated for the bag quarks j5^ 

The bag model has been used to calculate the GPDs in Ref . |49j , where the boosted bag 
wave function has been constructed using a simple prescription. In principle, one can perform 
a Fourier transformation of the GPDs to calculate the bag-model Wigner distribution. 

However, we choose a simpler way to calculate the Wigner distribution because the static 
bag has a fixed center. In fact, we can calculate directly from the wave function of quarks 
in the static nucleon just like in non-relativistic quantum mechanics. The rational for this is 
that assuming the GPDs are known, one can "correct" the relativistic effects associated with 
the boosted nucleon to obtain a Wigner distribution corresponding to the static structure, 
just like applying the relativistic corrections to extracting the static charge distributions. 
The Wigner distributions calculated from a static bag correspond to the ones with some 
relativistic corrections applied. 

As we have discussed in the last section, we define the Wigner distributions by the matrix 
elements of the Wigner operators W+(r, A;"*") in the hadron states. In general, because of 
translational invariance, only the off-diagonal matrix elements between nucleon states with 
finite momentum differences provide the 3D r dependence. However, in the static models 
such as the MIT bag, the quark wave functions are solved in the rest frame of the nucleon 
which has no translational invariance from start. With them, the Wigner distributions can 
be calculated as the diagonal matrix elements of the Wigner operator for the model nucleon 
fixed at the origin of the coordinates. For example, 

p+(f,x) = - / — e^^^(^ = 0|^(f-(A/2)n-)7+<^(f+(A/2)n-)|i? = 0) , (41) 

where \R = 0) represents the bag- model nucleon at R = and x = k'^/p'^ the light-cone 
momenta fraction of the proton carried by the quark, n a light-light vector with = 0,n~ = 
= 0. 

The quark field has the following expansion in the bag [5f 



^a{r,t) = N{nK){ba{nKm)tpnKj=i/2m{r,t) + dl{nKm)tp^n-Kj=i/2m{r,t)} , (42) 

n>0,K=±l,m 

where and dj^ are the quark and anti-quark creation operators in the bag, and N{nn) is 
a normalization factor. The wave function are solved from the Dirac equation with the bag 
boundary condition. For j = 1/2 and k = —1, one has 

^n,-i,i.(r, t) = ^( , ] e— . (43) 

For the lowest mode, we have n = 1, and cui^-i ~ 2.04. In the above wave function, a is the 
2x2 Pauli matrix, Xm is the Pauli spinor, and Rq is the bag radius, f represents the unit 
vector in the r direction, and ji are sphere Bessel functions. 

Substitute the above wave function into Eq. (jS)), we find the quark phase-space charge 
density, 

pi{r,x) = J ^e^^""'^^ [joiri)joir2) + ji{ri)ji{r2)fi ■ f2 

+i (jo(n)ji(r2)r^ - Jo(r2)ji(ri)ri)] , (44) 
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FIG. 3: The phase-space charge density p-\-{f, x) calculated in the bag model for values of Feynman 
momentum x = 0.1,0.33,0.5,0.9. 



where C/ is a flavor factor with Cu = 2 and Cd = 1 for up and down quarks, respectively. 
The position vectors fi and r2 are 

The above distribution satisfies the boundary constraint: Integrating over r yields the quark 
distribution function, while integrating over x gives the charge density of the quarks inside 
the nucleon. 

With p+(r, x), one can visualize the quark charge density as the function of x. In Fig. ^ 
we have shown a sequence of densities at x = 0.1,0.33,0.5 and 0.9. As the parton density 
indicates, the charge density is peaked around a; = 1/3 where the distribution is roughly 
spherical-symmetric. This is consistent with the finding that the bag model GPDs have a 
small ^ dependence. For smaller and larger x, the charge density can be negative. As x 
increases, the distribution at the center of the bag becomes smaller. As x further increases, 
the density there becomes negative. Similar phenomena happens as x decreases. Because 
the bag boundary limits the distance of the spatial correlation, the small-x distribution does 
not grow significantly as seen in experimental data. 



V. SUMMARY AND CONCLUSIONS 

In this paper, we have introduced the concept of the quantum phase-space distributions 
for the quarks and gluons in the nucleon. These distributions contain much more information 
than conventional observables. In particular, various reductions of the distribution lead to 
transverse-momentum dependent parton distributions and generalized parton distributions. 

Any knowledge on the GPDs can be immediately translated into the correlated coordinate 
and momentum distributions of partons. In particular, the GPDs can now be used to 
visualize the phase-space motion of the quarks, and hence allow studying the contribution 
of the quark orbital angular momentum to the spin of the nucleon. 

In light of this, measurements of GPDs and/or direct lattice QCD calculations of them 
will provide a fantastic window to the quark and gluon dynamics in the proton. 

We thank M. Burkardt and T. Cohen for many helpful discussions on the subject of the 
paper, and M. Diehl for critical comments. We thank the Department of Energy's Institute 
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by the U. S. Department of Energy via grant DE-FG02-93ER-40762. 
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